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A NEURAL MECHANISM FOR DISCRIMINATION: 
II. DISCRIMINATION OF WEIGHTS 


A. S. HOUSEHOLDER 
The University of Chicago 


A theoretical central mechanism for the discrimination of in- 
tensities as previously developed, together with plausible assump- 
tions concerning the receptors, are employed for the derivation of 
the discriminable difference between lifted weights as a function of 
the smaller of these weights. The function so derived depends 
upon three parameters, one parameter being the weight of the sup- 
porting member. Some empirical data are compared with the the- 
oretical predictions, and a few remarks are added to describe the 
physiological significance of the parameters. 


1. Some suggested central mechanisms for discrimination. In 
continuation of a previous discussion (Householder, 1939; this paper 
will be denoted hereafter by the letter D.) we consider a hypotheti- 
cal neural mechanism, whose elements have properties suggested 
by actual nerve fibers, and which has itself the property of diverting 
excitations of different intensities along different sets of paths. Our 
object will be to derive a formula for the just-noticeable difference in 
stimulus intensity for the particular sense of tension, and to compare 
this with some existing experimental data. 

For the most part, theorizing on the subject of discrimination of 
stimulus intensities is restricted, implicitly or explicitly, to the psy- 
chological level. Most satisfying of the purely psychological theoriz- 
ing in this direction are the psychophysical papers of L. L. Thurstone 
(1927, to cite only one), where the restriction is explicitly imposed 
and the argument unentangled by neurological concepts irrelevant to 
the problems and methods which he discusses. 

On the neurological level, E. D. Adrian (1932) notes the approxi- 
mately logarithmic relation between intensity of stimulation and fre- 
quency of response of a single receptor, an interesting point in con- 
nection with Fechner’s law. But evidently, on the one hand, a set 
of logarithmically responding receptors will not necessarily respond 
logarithmically as a set, and on the other hand, even if they do so the . 
nature of the discriminal process is still unsettled. 


t= 
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A. H. Holway and C. C. Pratt (1936) note the qualitative similar- 
ity between the curves for several modalities obtained by plotting the 
Weber-ratio against the stimulus intensity, and A. H. Holway and L. 
M. Hurvich (1937) conclude that “Differential sensitivity can be re- 
garded as a direct measure of the number of elements (= number of 
impulses per unit time) available for the discrimination of a minimal 
change in the intensity of stimulation.” In their discussion, Holway 
and Hurvich employ neurological terms chiefly to show the consis- 
tency of this notion with current neurological principles. No neural 
mechanism is described, however, to show how the available elements 
could be so organized that a lesser stimulus might traverse paths not 
reached by a stronger one. Oberlin’s (1934) suggestion that thres- 
holds of the receptor elements have a nearly normal distribution also 
falls short in this respect. Such an assumption, however, is to be in- 
cluded among those employed in the present paper (vide infra). 

H. Piéron (1933) makes a suggestion somewhat as follows. If 
two antagonistic effectors are reached by paths from the same centers 
with one having a higher threshold than the other, and if the one with 
higher threshold exerts some kind of an inhibiting effect upon the 
other, then indeed we should have a mechanism whereby quantitative- 
ly different stimulus intensities might lead to qualitatively different 
responses. Evidently the crucial features of this mechanism are the 
difference in thresholds at two synapses (or sets of synapses) and the 
inhibition exerted by one of these (sets of) synapses upon the other. 
Whether these elements be placed at the effectors themselves or at 
some intermediate station is of no importance theoretically. In a 
sense the mechanism here to be discussed is an elaboration of the one 
suggested by Piéron. 

N. Rashevsky (1938 a) and H. D. Landahl (1938) consider Sys- 
tems of neurons having the property that if two receptors (two sense 
organs or two parts of the same sense organ; e.g. two areas of the 
retina), O, and O,, are stimulated with unequal intensities S, and 
S,, and if stimulation of O, alone would lead to response f,, and 
stimulation of O, alone would lead to response R, , then the simulta- 
neous stimulations will lead to R, or R, alone according as S, exceeds 
or is exceeded by S.. If EF, and E, are the two effectors, it is only nec- 
essary to connect O, to HE, and also O, to E, by inhibitory neurons. 
With suitable quantitative relations the intensity of the response will 
depend upon the ratio S,/S, alone, and not upon the absolute magni- 
tudes of these intensities (Rashevsky, 1938 a, ch. xxvi). Landahl 
considers the same mechanism with the added assumption that either 
or both excitatory paths is subject to a random fluctuation. In a re- 
cently published paper Landahl (1940) has shown further that by 
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introducing suitable neuron circuits (Rashevsky, 19388 a, ch. xxiv) it 
is possible to delay the reaction to the first stimulus and to eliminate 
the requirement of simultaneous presentation of the two stimuli. 
Hence, knowing the distribution of the random fluctuations one may 
calculate the distribution of the judgments. If we call one of the 
stimuli the “standard”, presumably there will be a distribution curve 
for each standard, and with a suitable definition of a “noticeable” 
difference the Weber-ratio is defined for this particular standard. 

Another possible mechanism has been described by Landahl 
(1939), this one consisting of a set of neurons capable of transmit- 
ting excitations only if the stimulus lies in a limited intensity range. 

‘This consists of three neurons: an excitatory one and an inhibitory 
one, having the same threshold and leading from the same initial 
synapse, and both synapsing with the third neuron. The lower limit 
of the intensity-range is the common threshold of this initial neuron- 
pair; the upper limit depends upon the parameters of their excitation 
equations, and upon the threshold of the efferent neuron. Given a 
system of such neuron-triples with varying ranges, a stimulus S will 
act upon only those effectors reached by those triples whose range 
contains S. If, for example, we define two stimuli, S, and S2, as 
“completely discriminable” in case no neuron-triple of this system 
exists whose range includes both S, and S.,, then one could calculate 
the distribution of the ranges within the system from the empirical 
function defining, for the modality in question, the Weber-ratio in 
terms of the stimulus-intensity. 

The third mechanism, on which the ensuing discussion is based, 
is due, in its essential features, to Rashevsky (1938 b), and was de- 
veloped in greater detail in the paper D. Instead of accompanying 
each first-order excitatory neuron by a parallel inhibitory neuron, the 
inhibitory neurons interconnect the terminal synapses. Since the 
amount of inhibitory substance, 7, developed by the inhibitory neu- 
rons, depends upon the excess of « over 7 at these terminal synapses, 
and since this excess itself depends upon 7, we are led to an integral 
equation for the determination of « — 7 at each of the terminal syn- 
apses. If the primary neurons vary as to threshold it turns out, as 
with Landahl’s mechanism described above, that when any intensity 
S is applied to all the initial synapes, ¢ exceeds j at only a limited 
number of the terminal synapses. If two stimuli, S, and S., are de- 

fined to be “completely discriminable” when the sets of terminal syn- 
apses at which « exceeds 7 are distinct for these stimuli, then we can 

deduce the Weber-function defining the Weber-ratio in terms of the 


stimulus intensity. 
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To avoid a possible misunderstanding, let it be said that the def- 
inition of “completely discriminable” stimuli is not to be taken as im- 
plying that complete discriminability exists psychologically. Indeed, 
even at the neurological level J. Erlanger and H. S. Gasser (1937) 
find unpredictable fluctuations in the thresholds of individual neu- 
rons. It is enough to assume that stimuli which are just completely 
discriminable in our sense will be experimentally discriminated in a 
suitable percentage (say 75%) of the trials. 

Strictly speaking, another parameter should enter here into our 
Weber-function. Empirically 2 Weber-function could be obtained for 
each choice of the percentage required, and hence a one-parameter 
family of functions is empirically obtainable. It is assumed that one 
of these functions should agree with the function theoretically de- 
duced, but unless an additional assumption is made which will enable 
us to predict which:one it is (i.e. just what percentage corresponds 
to complete discriminability), this can only be decided experimentally. 
Actually we no not make this assumption, nor do we find experimental 
results sufficiently stable to allow the calculation of the family of 
curves with any reasonable accuracy. Hence in comparing theory 
with experiment we assume that the experimenter’s choice of the per- 
centage is correct within the limit of experimental error. 

2. Receptors and first-order neurons. We have spoken above of 
primary, or first-order neurons, initial and terminal synapses, and 
stimulus-intensities S. These terms were intended to refer only to 
the mechanism, and not to the nervous system as a whole. The “in- 
itial” synapses are initial for the mechanism, but are not necessarily 
peripheral. The stimulus-intensity S is the magnitude of « — j ap- 
plied to these synapses. The ‘terminal’ synapses are those at which 
the mechanism develops the « — 7 which excites the efferent neurons 
leading to the effectors. With this understanding, let S and S(1 + 6) 
be just completely discriminable. It is not necessary to repeat in this 
place the assumptions underlying the derivation of 6 as a function of 
S, since these can be found stated in the previous paper D. We shall 


merely assume the result contained in equations (D 15), (D 16), (D 
2a) and sDe23a): 


u=2,18/3, 
6=“+2, 
(1) 
ux + 7?—1=0, 


uzs—2(x+1)z— (4#+1)?=0. 
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The parameter 4 measures the “strength” of the inhibitory neurons 
acting among the terminal synapses. 

By assuming S to be proportional to the intensity of the peri- 
pheral stimulus, it was found that the theoretical curves agreed rea- 
sonably well with experimental data for light, sound, and touch. 
Moreover, only one parameter is involved, this parameter being the 
product of 4 and the factor of proportionality connecting S with the 
peripheral stimulus, R. It was pointed out, however (D, p. 52), that 
this assumption of proportionality would certainly require modifica- 
tion, and that some analysis of the peripheral mechanism would be 
required, in general, to yield the true Weber-function. The organs of 
. sight and hearing scarcely lend themselves to such treatment at pres- 
ent, for at least the intervention of visual purple, in the one case, and 
the organ of Corti, in the other, introduce unknown factors between 
the intensity of the adequate stimulus and the actual intensity of 
stimulation of the first-order neurons. The stimulation of the afferent 
fibers in muscle, however, seems to be much more direct, and in any 
case B. H. C. Matthews (1933) finds experimentally an approximate- 
ly logarithmic relation between the tension and the frequency of re- 
sponse. 

Suppose, then, that a weight is held motionless for a short time 
in, say, the open palm, the forearm being supported so that only the 
wrist muscles are employed in supporting the weight. The following 
assumptions are believed at least to offer no points of conflict with 
Matthews’ experimental findings: 

i) Each receptor has acting upon it a stimulus-intensity R 
which is proportional in magnitude to the weight being tested plus 

the weight of the supporting member, in this case the hand. 
; ii) The response of any receptor of threshold h leads to the de- 
velopment of S = « — j at the initial synapses of the discrimination 
mechanism according to the equation 


1 = log R/h. (2) 


iii) Thresholds of the receptors vary among themselves with 
a distribution function which is expressible either in the form 


N(h) =p h"(a—h)", (3) 
where n and m are positive integers, or, if the distribution is multi- 


~ modal, as a sum of such functions. Here we suppose the distribution 


unimodal. 
It follows that the total S developed at the initial synapses of the 


mechanism is 
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-R R R 
s= | S,N(h)dh =p | h(a —h)™ log — dh . 
0 0 


For convenience of calculation, note that a change of units for meas- 
uring R (and hence for measuring h) requires only a change in the 
value of the constants a and f’. Hence it is no restriction to take 
a = 1. This defines what we may term the “neurological unit” of 
weight, which is the weight which just excites the end-organ of great- 
est threshold. Its determination in a given series of observations 
amounts to the determination of one variable parameter. Thus, in 
place of the last equation we may write 


Rk R 
S=8' h(1—h)” log — dh. (4) 


“p 


If we expand by the binomial theorem and integrate we obtain 


Ty Sy OY Bee R? | 
SS PRS 2 (5 are 
and 
sap [wa —hym log dh 
; f log R se 
= iS as | Dp m 1 og 
pice C ec ta ete - 


By way of fixing the parameters m and n, it seems reasonable 
to suppose that the distribution function N(h) will be tangent to the 
h-axis at the limits of the distribution, i.e., at h = 0 and ath = 1. 
This corresponds to Oberlin’s assumption, referred to above, that this 
distribution should be nearly normal. We take the smallest integral 
values of m and n for which this is true, and therefore assume 


LR ees (7) 
which gives 
IN (hk) = Bh? kh)? (8) 
Then equation (5) becomes 
rere 
S = Te00 * (72k? — 225R -+ 200). (9) 


Since S is proportional to u by equation (1), then if we set 
“— 7 (10) 
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we have, after introducing a new parameter, 6 , 
12696? 7? = R? (72k? — 225Rk + 200) . (11) 


The reason for the factor 1269 will become apparent. 


For small 6 and hence large u we have, from (D386), with suf- 
ficient accuracy, 


2 
9 b] 
if we note that the z introduced in the above equation (10) is the re- 


ciprocal of that employed in (D36). Hence if A is the Weber-ratio 
corresponding to the peripheral stimulus R we have 


R*(1-+ A)* [72R2(1-+ A)? — 225R(1-+ A) + 200] 
= 1269 B* 79(1-+ 6) 


od =28451— (12) 


= 1269 P(e +224 21— 2) 
(13) 
2 
= 1269 6° (eas 


For R = 1 it follows from equation (11) that 27 6° 7? = 1. Hence 
let 
3pr=t (14) 
and obtain 
R?(72R? — 225R + 200) /47 = &, (15) 


R3(1-+ A)*[72R2(1-+ A) 2? — 225R (1+ A) + 200] /47 
= (t- 26)*—14 6°. 


We have assumed that 6 is small and u large in order to be able to 
write the approximate formula (12). Moreover, the last two equa- 
tions are only valid for R , R(1+4) and ¢ on the range from 0 tol. 
From equations (10) and (14) it follows therefore that 27/5? is equal 
to the reciprocal of the largest value of u on the range considered, and 
we may therefore safely neglect the term 146% on the right of the 
above equation. In fact in the experimental data to be considered later 
this term is of the order of 10-*, whereas the other terms are of the 
order of unity. Hence we have 


R3(1-+ A)*(72R2 (1+ A) ?— 225R (1+ A) 
+ 2001/47 = (t+ 28). - (16) 
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A still further, and final, simplification can be made as follows. 
Since equation (15) is valid only for R and t on the range from zero 
to unity we find that ¢ can be expressed as a function of R with suf- 
ficient accuracy by the interpolation parabola, 


t = 1.62R — .62R’, (17) 


which gives the correct values t = 0, .655, and 1 for R = 0, .5 and 1 
respectively. Similarly 


tt 26 =1.62R(1-+ A) — .62k?(1-+ 4)?. (18) 
Hence, on subtracting (17) from (18) we find 


28 = 1.62 R A — .62R? A(2-++- A), 
or 
31R? A? + (62k —81)RA+1006=0. (19) 
We choose the root of the quadratic (19) from the fact that 4d = 0, 
when £8 = 0, so that we have, finally, 


_ 81 — 62k — (81 — 62R)* — 12400 B 


2 
RA = : (20) 
or, if we set 
y = 124006, (21) 
81 — 62k — $1 — 62k)? — 
RAg ee ee eee (22) 


62 


Evidently if A is to be real for the entire range including R = 1 
we must have 


y < 861, B < 361/12400. (23) 


The minimum value of A is obtained from (19) by finding the small- 
est positive value of A for which the roots R of (19), expressed in 
terms of A, are real, i.e., the smallest positive A for which the dis- 
criminant of (19), regarded as a quadratic in R, is positive. This is 
24800 6 Iss 2 
6561 — 124008 6561 — » 
Evidently Anin increases in y(or 6). This minimum value is given by 
Bs 81 _ 6561— y 
opn62 (Ag eo yemerlo tae 


which decreases in y. Since the largest admissible value of y is 361, 


Ain = (24) 


Fin (25) 


A. S. HOUSEHOLDER 9 


the smallest value of Ruin is 50/81 = .62. For y = 0 we obtain the 
largest value of Rmin which is 81/124 — .65. At the same time Amin 
varies from its maximum value of 361/3100 = .12, for y = 361, to 
its minimum of 0 for y = 0. On differentiating the right member of 
(22) with respect to R it is easily seen that R A is monotonically in- 
creasing ink. 

i 3. The difference limen. If W is the weight of the stimulus-ob- 
ject in grams, W, the weight of the supporting member (the hand), 
then R is proportional to W + W,. We have then three variable para- 
meters, viz., W., B(or y), and this constant of proportionality. In 
complete form, we may write 


W+W,=uk, 
(W+W.)4 = 


81u — 62(W + W.) — v [8lu — 62(W + Wo)? — wy 
ee oa (26) 


where )W + W,)A is the experimentally measurable difference limen 
AW. However, W, is also experimentally measurable so that only 
two of the parameters, u and y, are incapable of direct determination. 

Nevertheless W, is not ordinarily given, although one can cer- 
tainly tell its order of magnitude from the weight of one’s own hand. 
We present herewith a comparison of the theoretical curves with some 
experimental data obtained by A. H. Holway, J anet E. Smith and M. 
Zigler (1937) for unimanual weight discrimination. No attempt has 
been made to fit the data given by Oberlin (1934), which seem rather 
less regular, although Oberlin’s curves seem to follow the same gen- 
eral course. Unfortunately in their experiments, Holway, Smith, and 
Zigler seem to have made no effort to rule out possible pressure cues, 
so that actually more than muscular sense is involved. Data are de- 
termined for the two hands separately by the three observers, giving 
six sets in all, for weights ranging from 100 g. to 1600 g. We may 
note that in graphing their own data these writers take the mean for 
the two hands of each observer, thus obtaining much greater regu- 
larity than is to be found in the separate sets. | 

We present, in the accompanying figure, the graphs of the data 
for the separate hands of each observer, as well as the graph of the 
mean for each observer, along with the graphs of the best fitting theo- 
retical curves. The reason for presenting the graphs of the means is 
the following. If the sensitivity of the two hands is essentially inde- 
pendent of the “handedness” of the individual, the data for the two 
hands should be equivalent and the means would thus be more reliable 
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since they were obtained from twice as many observations as were 
the data for the hands separately. And indeed, in the cases of ob- 
_ servers 1 and 2 the greater regularity is at once apparent. We note 

further that the means for these observers also fit the theoretical 
curves much more closely than do the data for the separate hands. 
This is not true, however, for observer 3, where the points are equally 
irregular and the fits equally uncertain for the means and for the 
separate hands alike. 

A mention of some of the steps in the determination of the para- 
meters is relevant. Since the hand weights were not given, W. was 
determined along with uw and y from the data for each hand separately, 
so as to obtain the best fitting theoretical curve. The theoretical curve 
was fitted to each of the six sets of data independently of the others, 
so that when the values of the parameters had been estimated for 
one hand of one observer, no account of these values was taken when 
the parameters for the other hand of the same observer were being 
estimated. Nevertheless, within the limits of accuracy allowed by the 
data, the values of W, for the two hands of the same observer came 
out the same in every case. Moreover, it turned out that for both 
of observers 1 and 3, the estimated value of W, was 400 g., whereas 
for observer 2, W. came out to be 350 g. Since observer 2 was a wom- 
an and the other two men it seems likely that she should, in fact, have 
smaller hands. 

In fitting the curves to the averages, the values of W, were taken 
as known from the previous determinations, but « and y were deter- 
mined independently. With each graph is given the values of the three 
parameters. The label O,, means that this is the graph for the left 
hand of observer 1; Ou.» signifies the mean for the two hands of ob- 
server 1; and similarly for the others. 

It was pointed out above that the parameters uw and y are not sub- 
ject to a direct determination. However a word might be added by 
way of pointing out their physiological significance. It will be recalled 
that R is the ratio of the applied stimulus (hand-weight plus object- 
_ weight) to the maximum threshold of the sensitive elements. For 
R=1,W+W.=#,s80 that « — W, is the least weight necessary to 
excite all the afferent fibers in the group utilized for the discrimina- 
tion. It may be, however, that another group of fibers is present, pos- 
sibly itself distributed according to a function of the form (3) but 
with different parameters and having a mode greater than the mode 
for the first group — greater, in fact, than » — W,. If so the 4W 
curve, which is rising rapidly where W=n—W., will level off some- 
what for larger values of W. If not, either the increase in W will con- 
tinue, and at an increasing rate, or else in effect the judgment will no 
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longer be made with the wrist muscles alone and additional cues will 
be employed by the subject. These are questions for the experimenter 
to answer. 

Whatever the case may be, the foregoing theoretical considera- 
tions indicate the need for the determination of additional difference 
limens under the following conditions: 

1. When the wrist and finger muscles are not acting, nor yet the 
shoulder muscles, but only the elbow muscles. 

2. When only the shoulder muscles are acting. In such experi- 
ments a plaster-of-paris cast such as the one employed by M. von Frey 
(1913) would doubtless be needed, so that purely tactile cues could be 
completely eliminated. 

As for 6 (or y), this is inversely proportional to the cube root of 
the value of wu, equation (1), corresponding to R = 1. Hence the 
smaller the value of 8 the smaller the value of 6 at R = 1, and thus 6 
gives a measure of the fineness of the discriminal mechanism, a small- 
er § meaning a greater precision. 
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BIOLOGICAL SELF-REGULATION AND OF CELL POLARITY 
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In a cell, in which the permeability to a metabolite is a function 
of the concentration of that metabolite, situations may occur, in 
which the diffusion field will exhibit certain assymetric patterns, 
even though the cell may possess geometrically spherical symmetry. 
This pattern results in a polarity of the cell. Moreover, the pattern 
being the result of a dynamic equilibrium, it possesses the property 
of self-regulation. Dividing the cell in two results in the appearance 
of a similar pattern in each half-cell. 

Another case when such self regulation and polarity occur is 
given by considerations of the action of the diffusion forces upon 
colloidal particles, which affect catalytically the metabolic reactions. 
A simple case is treated mathematically. 


The problem of biological self-regulation offers many different 
aspects. One of the most fundamental of these, and at the same time 
one of the most difficult is that of the development of regulation eggs. 
Most drastic interference with the visible internal structure of these 
eggs may leave them still capable of producing normal embryos (Har- 
vey, E. B., 1932, 1936). Such self-regulation has perhaps been the 
phenomenon, upon which claims of the existence of non-physical ele- 
ments in biology have been based. While much has been written on 
the physical aspects of self-regulation and it has been pointed out, 
that even relatively simple physico-chemical systems may at times 
possess the property of self-regulation, yet all these general consid- 
erations and analogies do not bring us much nearer to the solution of 
the problem. If we are to incorporate eventually these phenomena 
into the field of mathematical biophysics, we must proceed by study- 
ing theoretical models, which are sufficiently specific to be applied not 
only to the phenomenon of self-regulation in general, but to some def- 
inite biological case of it. Of course, when such a model is developed, 
it is by no means to be regarded as the explanation of the particular 
phenomenon, for which it-is intended as a model. Most likely there is 
a large number of different possible explanations. It is only by study- 


ing systematically all theoretically possible cases and by a comparison 
sith RS 
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of the mathematical deductions from each of those cases with avail- 
able experimental data, that we may in the future arrive at a decision 
as to which, if any, of the theoretically studied situations actually 
occurs. 

With this in mind, we shall now investigate a problem, which is 
suggested by such observations, as for instance those by H. B. Harvey 
(1932, 1936) on centrifuging Arbacia eggs. Subjecting them to a 
force 10,000 times that of gravity not only produces a definite strati- 
fication of the content of the egg, according to the specific gravity of 
the different constituents, but results eventually in a breaking up of 
the egg in two and even in four parts, each part containing different 
visible components. Yet, when fertilized, all four parts can give rise 
to normal embryos. 

The early differentiation of the fertilized egg indicates definitely 
the preexistence of some sort of nonuniformities. At the same time 
these nonuniformities obviously cannot be connected with or even re- 
lated to the visible material. Moreover, whatever latent nonuniform- 
ity pattern may have existed in the normal egg, the breaking of the 
latter into four parts of very different gross constitution must have 
resulted in each part acquiring again the same latent pattern, which 
characterized the original normal egg. 

A suggestion has been made that the pattern responsible for the 
early differentiation is contained not in the cell volume, but on the 
surface in the cell membrane (Weiss, 1939). The latter being of a 
much more rigid nature than the more liquid interior of the egg, the 
centrifugal forces used in the experiment may not be sufficient to up- 
set the inner structure of the membrane. Since there is a direct mi- 
croscopic evidence for a superficial localization of the early noticeable 
differentiation of the egg (Weiss, 1939), it appears of interest to 
investigate this possibility more quantitatively. 

It must however be remembered, that while there is good evi- 
dence for a more solid-like consistency of the egg membrane as com- 
pared with the interior of the egg, yet the membrane certainly cannot 
be considered as a solid bag, containing a liquid interior. For in such 
a case a force strong enough to break the egg into parts would simply 
tear the solid membrane, leaving the interior to flow out. The micro- 
photographs taken by E. B. Harvey (1932, 1936) and reproduced by 
P. Weiss (Weiss, 1939, p. 191) show that under the influence of the 
centrifugal force the egg elongated and then divided by constriction 
in the middle. In other words the membrane at the equatorial region 
fuses together as only a liquid would do. If however the membrane 
is merely a very viscous liquid, sufficiently strong forces would upset 
its internal structure also. 


N. RASHEVSKY 17 


I 

Let us consider a spherical cell of radius 7, cm., producing a sub- 
stance at a constant rate gq gm cm™ sec. Let the external concentra- 
tion of the substance be c, gm cm~ and the external diffusion coef- 
ficient D, be very large, so that we may put D, = o. This merely sim- 
plifies the calculations, without introducing any essential limitations. 
Let the internal diffusion coefficient be D cm? sec*. Assume further, 
that the membrane of the cell is not uniform in its physical constitu- 
tion and that therefore its permeability is not the same at every 
point. Let the structure of the membrane possess an axial symmetry, 
thus imparting a polarity to the cell. Let at one pole the permeability 
have the value h, cm sec”, while at the other pole let that value be 
h. em sec. Due to such a nonuniform and asymmetric distribution 
of the permeability, the distribution of concentration of the produced 
substance within the cell will also not be spherically symmetric. In 
particular the concentration inside the cell at the membrane will vary 

from point to point. 
To calculate the average concentration distribution in such a cell 
we shall use the usual approximation method (Rashevsky, 1940). Let 


h, denote the average permeability in one hemisphere (Fig. 1), while 


FIGURE 1 


h, denotes the average permeability in the other one. Correspond- 
ingly we shall denote by ¢: and c, the average concentrations at the 
membrane inside the cell in each hemisphere. Denoting by ¢ the aver- 
age concentration within the cell, we have 


2D <= In (ox — 0) § (1) 
0 
2D <—* = hae —%) § (2) 
0 
— e—c 
5are'd = 2Qar? X 2D ~ . 4 Onn? X 2D =a (3) 
j 0 9 
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Equations (1) and (2) give: { 
2D¢ =e, 2D6 ASrehse, 
eee ee ome 0. es eee 


Cita ES ? 2 rs (4) 
2D + 7h, 2D + rhe 
Hence ated. 
h.(¢ = —C¢ 
ple epee ACCC ated 5 sige Tol NSgetalyy (5) 
2D er ha 2D + rohy 


Introducing equation (5) into equation (3), we find after rearrange- 
ments: 


To (2D rohx) (2D + Toho) 
3. Dh 2DS_ Ri yee hs (2D Beh ye 


C= e+ (6) 


For h, = ho =h , that is for a uniform membrane, equation (6) re- 
duces to 


6D’ (7) 


as should be the case*. 
Equation (6) introduced into equation (4) gives 


2 2D + ols : 
Cy — Co + ee, SO 9 
3 h,(2D4- rhe) +. hy (2D + rohy) 
‘ (8) 
2D -L rah 
os = 0 + 2rd _ 2D +40, 


3 hy (2D 4 rhe) + hy (2D 4 rh) 


Thus an asymmetry of the cell membrane will result in an asym- 
metry of the diffusion field within the cell. For the approximate cal- 


culation of this field we considered two discrete average values h, and 


h. of the permeability. The actual value of the permeability is as-_ 
sumed to vary continuously along the meridian. If this variation is 
monotonic, then a division of the cell into several] parts by sections 
parallel to the equator will produce smaller cells, each having a similar 
polarity. 

We have here a model which may reproduce some of the prop- 
erties of actual cells. Yet this model does really. not show any regula- 
tion. It is essentially a “mosaic” model, the mosaic being attributed 
to the membrane and possessing a rather simple pattern. The above 


* The equations derived for a sphericall 
evsky, 1940) differ from equation (7) in the values of the numerical: coefficients, 
which are respectively 2/9 and 1/9. This is due to the fact that we considered 
there actually a cylinder whose length equals its diameter. 


y symmetrical cell previously (Rash- 
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considerations suggest however an interesting possibility of actual 
regulation, which we shall mention here. 

Suppose that the permeability h of the membrane is itself in- 
fluenced by the concentration ¢ of the metabolite in its vicinity. Such 
a supposition is biologically rather plausible. Let this influence be 
such, that an increase of c results in a decrease of h. Let the perme- 
ability be originally uniform all over the membrane. Then the diffu- 
sion field will possess spherical symmetry. Let now the concentration 
become, because of some disturbance, slightly higher at one region of 
the membrane. This will result is a decrease of permeability in that 
region, thus creating a condition similar to the one discussed above. 
But in that case the concentration c at the membrane will be larger 
in the region with smaller h. This can be readily seen from equations 
(8), which give: 


27.2q (le — hs) 


Cpt ee (9) 
hy (2D + rola) + hy (2D + Toh) 


—— 


Equation (9) shows, that when h, < ho, then c, > c.. Thus the re- 
sulting change in permeability will result in a further increase of ¢ at 
that place, thus enhancing the asymmetry. We may ask whether un- 
der such conditions stable asymmetric configurations may not be ob- 
tained. 

Mathematically the problem reduces to the following. Let h bea 
prescribed function of c, h = f(c). We then have h, = f(¢,) and 
h. = f(¢.). Denote the right hand side of the first equation (8) by 
u,(hy, he) and the right hand side of the second equation (8) by 
U» (NR, he). Then, since h, and h, are functions of ¢ and c., we find 
that u, (Tu, ie) = ulf (cr), f (C2) ] = V1 (Gr, C2) ;and similarly w(t, he) 
= ,(C,, G2). We must now look for stable solutions of the system of 


equations: 


Cy = V1 (C1, Ce), 
(10) 


C3) = V2 (C1, C2), 


other than c, = ¢ , which corresponds to a spherical symmetry. It can 
be shown that when f(c) is monotonic, only spherically symmetric so- 
lutions are possible. Physically this is due to the fact that when an 
asymmetric configuration is obtained in a manner described above, 
then any slight increase of c at the lower end will result in a decrease 
of h at that end, with a resulting further increase of c, until it reaches 
the same value as at the higher end. However we cannot @ priori ex- 
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clude the possibility that for some special types of the function f(¢) 
asymmetric solutions of the system (10) may exist. The problem now 
is to investigate whether this is possible and if so, for what forms of 
f(c). Furthermore we shall have to ascertain, that if such possibilities 
exist, the necessary form of f(c) is biologically plausible or at least 
possible. iy 
For given h, and hy , the values ¢, = u, (Ai, ho) and C2 = U2 (My, he), 
as expressed by equations (8), are obtained by considerations of ma- 
terial balance in a stationary state. Therefore when ¢, > U,, ¢, will 
decrease until it reaches the value u,. Conversely if ¢, < u,, ¢, will 
increase until it reaches the value u,. A similar thing holds for ¢. 
Therefore the conditions of stability require that any small increment 


of c, would result in a lesser increment of u,(h,, he) or, which is the 
same, of ¥,(¢,, ¢2), so that dv,/dc, < 1; and similarly 0v./oc,. << 1. 
For the symmetrical case, when c,; = c,, it may happen that the con- 
figuration will be stable for a simultaneous equal increase or decrease 
of c, and ¢,, but unstable for an increase of c, accompanied by a de- 
crease of ¢, or vice versa. 

If asymmetric stable solutions of the system (10) exist for some 
forms of f(¢), and if they exist for a sufficiently large range of values 
of 7), then we have a case of a true self-regulation. No matter how 
we divide such a cell, whether equatorially or meridionally, no matter 
how we stir up its content, the cell or each part of it within the ad- 
missible range of 7,, will automatically reconstitute the asymmetry 
of the diffusion field. 

We have however rather formulated a problem than given a so- 
lution. A consideration of a different physical mechanisms leads us 
to a simpler solution of this type of problem. 


II 


Consider again a spherical cell, producing a substance. Let how- 
ever now the permeability be a constant, h . Let moreover the cell con- 
tain a catalyst, which inhibits the reaction, so that the higher the con- 
centration n of the catalyst, the smaller the rate of production q. If 
the molecules of the catalyst are rather large, or if the latter is pres- 
ent in a colloidal state, with a particle size > 2 X 10-7 em., then as we 
have seen elsewhere (Rashevsky, 1938), due to the action of the diffu- 
sion forces, these particles will show a pronounced nonuniformity of 
distribution, their concentration n being larger in regions of smaller 
concentrations ¢ of the diffusing metabolite. For a produced sub- 
stance, when the negative gradients are directed outwards, the cata- 
lyst will accumulate at the periphery, decreasing there the rate of pro-. 


+ eA 3 ~ 
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duction and still enhancing the gradient of c. However, for a spheric- 
ally symmetric distribution of the catalyst, everything will remain 
spherically symmetric and the average rate of production q will be 
the same in any arbitrarily chosen half of the cell. If however, for 
some reason, in one of the hemispheres the average concentration c 
increases, this will result in an increase of n in the other hemisphere 
at the expense of the first. This will cause a further increase of ¢ in 
the first hemisphere, enhancing the asymmetry. To investigate the 
possibilities in this case, let us first investigate what effect a variation 
of gq from point to point in a cell will have upon the diffusion field. 

Referring again to Fig. 1, denote the average concentration in 
one hemisphere by c,, in the other by c.. Denote the average rate of 
production in the first hemisphere by q, , in the second by q.. We now 
have 


eT Cy 5) ; (11) 
Ne 


Cy — Co 


2D = (ts —"C5) « (12) 


Lo 

The material balance equation is now however somewhat modified. 
Due to the different values of c in the two hemispheres, there will now 
be not only an outflow of the metabolite from the cell, but also an in- 
ner flow, from one hemisphere to another. Since the average gradient 
for this flow is approximately (¢c, — ¢:)/7ro, and the area through 
which the flow takes place is 27,2, that of the equatorial circle, the to- 
tal inner flow is equal to 27,D(c,—c¢.). We now have for the first 
hemisphere 


ae ap a! is 
2 seheqy = Bary? X BDA + ar D(C —) , (43) 
0 
and for the second 
== 6, = = 
= a = 2ar? X 2D as =~ A od) (07-63) ss (14) 
0 


Equations (11) and (12) give 


2De,+ Toheo ees 2Déo + TohCo 
Saipan y iter Be) 2b. rh 


Introducing (15) into equations (13) and (14) we find after rear- 
rangement: 


C (15) 


22 MATHEMATICAL BIOPHYSICS 


3D (2D + 5r,h)¢, = 83D (2D + rh) 
+ 12Dr he) + 2(2D + 1h) 107: 5 


(16) 
3D (2D + Breh) Gy = 8D (2D + Toh) &, 
+ 12Dr,he, + 2(2D + rh) 1o7Qr « 
Solving we find: 
ie i 19° (2D + rh) [ (2D + rh) de + (2D + 57r,h) q;] 
eanitis= 12Dr,h (2D 3r.h) , 
(17) 
— _, , rat(2D + rh) [2D + 7eh) dy + (2D + 5rd) 
pace 12Dr,h(2D 3r.) ats 
Hence 
ee i 0” (2D rh jeer e 
sales Spe Tol (Be teat Ces ta) (18) 


83D (2D + 38r,h) 
The difference c, — ¢, vanishes for gq; = q2. When q, > q,, then 
CE yee. 

Now let us investigate the effect of such an asymmetric concen- 
tration distribution of the metabolite upon the concentration n of the 
catalyst. Denoting by V the volume of the particle, by N — Avogra- 
dro’s number, and by M — the molecular weight of the diffusing meta- 
bolite, and putting 


(19) 
we have (MB, p. 67), for the ratio n,/n. of the concentrations of the 
catalyst in the two hemispheres 


a = e-4(61 - 2) : (20) 


Denoting by ” the average concentration of the catalyst in the whole 
cell, and remembering that the total amount of the catalyst 2 arn is 


constant, we have 


(21) 


Putting 
(%—C,=72, (22) 
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we find from equations (20) and (21) 
Ns. —N, = 2n tanh fax. (23) 
For the effect of the catalyst on the reaction rate q let us con- 
sider the simplest possible relation, a linear one. 
Let 
q=%—an, (24) 
where a is a constant. Then 
G1 — G2 = (M2 — M1) « (25) 


If an asymmetric distribution is to be possible, it must be represented 
by a stable root of the equation, obtained by introducing (23) into 
(25) and then introducing the latter into (18). Putting 


To2(2D + rh) ¥. 


SDLP ore eS) 
we thus find 
x = 2Aan tanh fax. (27) 


If, as will usually be the case, ax is rather small, we may expand the 
hyperbolic tangent keeping only the lowest nonvanishing power. Equa- 
tion (27) then becomes 


1 
= Se aa ogre 2 
¢ = Aanax (1 3°) , (28) 


which, besides x = 0, has a positive root: 


nie 3 (Aana — 1) 
provided 
Aana>1. (30) 


That the root (29) corresponds to a stable configuration is seen by 
graphing equation (28). The curve which represents the right hand 
side of (28), is nothing else but the right hand side of equation (18). 
By a similar argument as before we must have for stability 


d 1 22 = * 
rae 1 ft Ti og 
qq Aonon (1 70) < or 


which is actually the case. 
The model discussed here does exhibit true selfregulation. No 
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matter how the cell is divided or how its content is stirred, the char- 
acteristic asymmetry of the diffusion gradient will be automatically 
reestablished. By considering the presence of several catalysts with 
different particle size, and affecting different reactions, we may ob- 
tain complex asymmetric patterns of the diffusion field, which thus 
may represent a rather intricate latent self-regulating structure. 

It must be noted that the diffusion force per unit volume, 
—(RT/M) grad c may be of the order of 107 — 10% dyn. cm’, (Rash- 
evsky, 1938), which is not only comparable, but even may exceed the 
centrifugal forces used in the above mentioned experiments. More- 
over while the effect of the centrifugal forces upon a particle is pro- 
portional to the difference in specific weights of the particle and of 
the surrounding liquid and vanishes when the two specific weights 
become identical, no such restriction holds for the effect of diffusion 
forces. The model discussed here may in principle be mutilated be- 
yond the point when self-regulation is possible, by applying such a 
strong centrifugal force, that all the catalyst is driven into one cell 
fragment. But this may require centrifugal forces far in excess of 
those available now. 


III 

Another important result of this study is that it leads us into the 
theory of cell polarity. A cell like the one discussed in section II al- 
ways possesses a “polarity,” and it is known that polarity is rather 
the rule in cell biology. The asymmetric distribution of the metabo- 
lite and therefore of the diffusion forces will in general result in 
asymmetries of all other important physical properties of the cell. 
Aggregates of such polar cells will themselves exhibit polar proper- 
ties and a way is thus indicated for a mathematical biophysics of dif- 
ferent types of tissues. 

A polarity considered in section II is of a dynamic nature. But 
due to the irreversibility of some biological reactions, it may result 
after a lapse of time in a permanent polarity of a static nature, as 
discussed in section I. 

The author is indebted to Dr. A. S. Householder for a discussion 
of this paper and checking of the calculations. 

The publication of this paper was made possible by a grant from 
the Rockefeller Foundation to the University of Chicago. 
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; The action of diffusion forces in aggregates of metabolising cells 
is studied mathematically, and it is shown that under definite con- 
ditions these forces may lead to the formation of inner cavities 
inside of the cellular aggregate. Different quantitative relations are 
derived and the possible bearing of these results on some embryologi- 
cal phenomena is discussed. A possible application of these consid- 
erations to the theory of formation of vacuoles is also discussed. 


In our previous studies (for a summary see Rashevsky 1938 and 
1940; hereinafter referred to as MB and A, correspondingly) we have 
discussed various dynamic aspects of the action of metabolic diffusion 
forces upon the cell. We have only occasionally touched the dynamics 
or even statics of cellular aggregates. In this paper we propose to 
study some of those omitted questions. 

We have seen that in any cell with positive effective metabolism 
(A, chapter iii), the diffusion forces are directed outward and in the 
case of any deviation of the cell from spherical shape, a situation 
which is always present, those forces result in the elongation and 
eventual division of the cell. As has been remarked in MB, chapter 
viii, in principle we may consider also those effects of the forces which 
tend to simply disrupt, to burst, or to expand the cell. If we consider 
the cell as we have done hitherto, as a-liquid system, then due to the 
practical incompressibility of liquids such an expansion would re- 
quire forces much larger than those which may be produced by dif- 
fusion. If on the other hand we consider the expansion of the cell as 
merely the expansion of its membrane, followed by a free flow of 
water from outside into the increased volume bounded by the mem- 
brane, then such an expansion becomes quite possible under the ac- 
tion of diffusion forces, and its study leads to some interesting as- 
pects of cell growth (A, chapter iv). 

A somewhat different situation is possible when we consider not : 
a single cell, but an aggregate of cells. If such an aggregate has a 
positive effective metabolism, then the individual cells of the aggre-. 


A, ae 
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gate repel each other. If nothing holds these cells together, then they 
simply would separate and disperse under the action of the diffusion 
forces. If however the aggregate is held together, for instance by an 
elastic membrane in which it is enclosed, or by forces of surface ad- 
hesion, which makes the cells stick to each other, then this membrane 
will oppose a complete scattering of the cells. The cells may however 
move apart to some extent, stretching the membrane. Since the vol- 
ume of each individual cell remains constant, such a moving apart 
will result in the formation of a cavity inside the aggregate. If the 
aggregate is approximately spherical, this will result in the formation 
of a hollow spherical shell. Such a shell need not consist of a single 
layer of cells. The thickness of the shell and its overall size will be 
determined by the equilibrium between the diffusion forces of repul- 
sion and the elastic forces of the membrane, which tend to counteract 
the forces of repulsion. 

To establish the conditions of equilibrium, we must first of all 
calculate the work done by the diffusion forces when the shell ex- 
pands, becoming thinner. Denoting by 7, the radius of the cavity, by 
7’, the external radius of the shell ,we see that in such an expansion 
both 7, and 7, increase. Since however the total volume V of all the 
cells remains constant, we have 


4 
gts ae? =V, (1) 
Putting 
BV; 
OREN Se: (2) 
we may write 
To? —=Ts° —, (3) 
or 
Te N) Nar ily (4) 


On each element of volume of the shell there acts a force (A, 
chapter iii) 


8 RT 
F=— —-——-gradc, 
2M : (5) 


the notations being the same as in A. Moreover, at the outer surface 
of the shell there acts a force due to the difference of concentrations 
on both sides of the membrane. Without introducing any fundamen- 
tal restrictions, we may consider a case in which this surface force is 
negligible, due to a large permeability of the membrane. This merely 
simplifies the final expressions. 
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The total work done in any deformation of the aggregate by the 
force (4) is equal to (MB, p. 74): 
= 3 RT u Vv 
OW = Se= ) Vk : 0) grad c,dt ’ (6) 
where b is the displacement of the element of volume dz and c; the in- 


ternal concentration. But 
ii 


Sf vgrad cdr = — f, e,divodr+ f, csv dS , (7) 


where c, denotes the concentration at the surface which bounds the 


volume V of the aggregate, » denotes the value of » at the surface and 
dS is the element of the surface. 

Since in the deformation of the type considered here not only the 
total volume V remains constant, but also the volume of every indi- 
vidual cell is considered as constant, div » = 0, and (5) and (7) give: 


Leelee ier. ies eae Cf) 

A spherical shell is bounded by two surfaces: the inner sphere 
of radius 7, and the outer one of radius 7,. The integral in (8) con- 
sists therefore of two parts. Denoting the concentration at the inner 
surface by c, , that at the outer by ¢,, and remembering that the out- 
ward directions of the normals at the two surfaces are opposite and 


that therefore at the inner surface » = — dr, , while at the outer sur- 
face b = dr., we find that the integral in (8) is equal to 
— Anr,?c,dr, + 4nr,?cod7rz . (9) 


Since during the deformation considered the volume of the shell re- 
mains constant, 4ar,2dr, = 4zr.?dr.. Hence, combining (9) with (8) 
we obtain: 


62kT u 
OW = a 122 (C1 — C2) dT. . (10) 


The quantity (c, — c.) may be easily obtained for a perfectly spheri- 
cal shell by integrating the diffusion equation. Since however there 
are no perfect spheres anyway, we shall use the usual approximation 
method (A, chapter i) and save considerable calculations. 

Since we consider the case of positive effective metabolism, we 
have c, > ¢, and the total ouflow is equal to 


Anr.?D; (¢, — C2) 
pts, 


(11) 


? 
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D; denoting the internal diffusion coefficient. 
The total production equals for a case of constant rate of produc- 
tion: 


Vg= snag. (12) 


Putting (11) equal to (12), solving with respect to c, — ¢, and sub- 
stituting into (10) we have 


is 2aRT wga (12 — 11) 
MD; 


Now we must calculate the work done by the elastic forces of the 
membrane, when it is stretched due to an increase of the radius of the 
shell from 7, to 7. + dr.. Assuming the membrane to behave like an 
elastic body, the work done will be equal to edS/S , where ¢ is a con- 
stant similar in nature to the modulus of elasticity. Since S = 4a7”,?, 
dS = 8ar.dr, therefore the work done by the elastic forces of the mem- 
brane is equal to 


OW dr2. (13) 


Wipe oe (14) 
Ye 
The shell will be mechanically in equilibrium, when the sum of (13) 
and (14) is zero. Writing down this requirement, substituting in (13) 
for 7, its expression (4) and putting 


Re at (15) 
we find 


12 (2 — Wri—a)=K. (16) 


For any given value of a, or what amounts to the same, for any given 
value of V , equation (16) gives us the corresponding value of 7. and 
because of (4), also the value of 7,. Unfortunately equation (16) is 
not conveniently solved in closed form. A discussion of it is however 
of interest. Physically 7, cannot be smaller than ¥/a. This corre- 
sponds to 7, = 0, that is to the absence of any cavity. If for a value 
of a the largest root of equation (16) is less than \/a@, this means that 
for this value of a the stable configuration is that of a compact aggre- 
gate of cells, without any other inner cavities. This is due to the fact 
that in this case the volume of the aggregate is still so small, that the 
gradients of concentration and therefore the diffusion forces are not 
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large enough to overcome the elastic forces of the membrane. As the 
volume of the aggregate increases, these forces increase, and for a 
critical volume V* , to which corresponds the value a* of a, the diffu- 
sion forces will just balance the elastic forces. For this volume we 
have r.> = @, or 7, = Wa. As the volume increases further, a cavity 
is formed andr, > Wea. 

The value a* can be found by putting in equation (16) atau wa 
and solving the latter for a. We then find: 


5 M Die 
a* = K3/5 = ( = 3k 17 
RT Gn ) (17) 
The outside radius 7» is in this case equal to 
pa ead ade (18) 
RT qu 


Equation (18) gives us the size of the cell aggregate above which it 
will not remain closely packed, but will form a hollow shell. 

Instead of varying the size of the aggregate we may keep it con- 
stant and vary any one of the other parameters, for instance q. AS 
long as 72 < WK, the cellular aggregate is a closely packed one, there 
are no cavities and 7, = 0. If we keep 72 = 7 constant, but increase 
q, then K will decrease and eventually we shall have 7 > WK. As 
soon as this happens a cavity will be formed, and 7 = 7 will increase 
due to the expansion of the whole structure. The size of the cavity 
will be determined by (4) into which the positive root of equation 
(16) is to be substituted for 7. 

If instead of an elastic membrane we consider the outside mem- 
brane as a liquid film with the surface tension y erg cm™ we obtain a 
somewhat different result. The work done by the surface tension is 
equal to — y dS = — 8ayr.dr,. Using this instead of (14) we now 
find instead of (16): 


- (19) 
RT qu 


a(t,—Wre—a) =K'r, ; K’ 
Equation (19) may be written 
(a— K')r, = aw ro — a. (20) 


Raising both sides to the third power and solving for 72, we find: 


ae? a* 
To =a . (21) 
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Since again 72° SS a, we find the critical value a* for which 7.3 > @ 
from 


at aon 
which gives 
Ga hes (22) 
The corresponding value of 7,* is 
714M D, 
r* = WK = ee sy. (23) 
RT qu 


an expression analogous to that for the critical radius of a cell, above 
which it divides (MB, chapter viii; A, chapter iii). 


II 


As another interesting theoretical possibility we shall investi- 
gate a case of formation-of cavities due to the inter-play of diffusion 
forces only, in the absence of any elastic membranes. 

Consider again a densely packed cellular aggregate, consuming 
at a constant rate per unit volume oxygen supplied from outside. The 
concentration of oxygen decreases from periphery to the center and 
cells closer to the center have a smaller “external’’ oxygen concentra- 
tion, than cells closer to periphery. 

It has been shown by H. D. Landahl (1939) that the glycolytic 
coefficient of a cell (A, p. 34) increases with decreasing external 
oxygen concentration (A, p. 36). On the other hand, we have shown 
that the effective metabolic rate is a linear function of the glycolytic 
coefficient 6, of the form 


q= Ap "By A> 0. Sah 2 Oe (24) 
When B < 0, q is zero for a positive value of 8 = — B/A, positive 
for 6 > — B/A and negative for 8 < — B/A. Buta change in the 


sign of gq means a reversal of the direction of osmotic forces, which 
are forces of attraction for g < 0 and forces of repulsion for q > 0. 
(A, chapter iii). 

With a sufficiently large supply of oxygen to the cellular aggre- 
gate, the cells nearer to the periphery may have a sufficiently small 8 
to give gq < 0, and will attract each other, while the inner cells, with 
a smaller supply of oxygen, may have a large enough £ to make q > 0, 
and will repel each other. Thus within the same aggregate we shall 
have a system of two opposite forces, one tending to expand the agegre- 
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gate, the other tending to keep it together. In a sense we may Say, 
that the outer layer of the mutually attracting cells now plays the role 
of an elastic membrane, that holds together the inner cells, which re- 
pel each other. If the repelling forces prevail the aggregate will ex- 
pand, forming a cavity. But due to this expansion the thickness of the 
layer of cells, forming the hollow shell, will decrease. Therefore the 
inner cells will receive more oxygen than before, their glycolytic co- 
efficient 8 will decrease and the repelling force will decrease. An equi- 
librium is reached, when the two forces just compensate each other. 

To determine the equilibrium configuration we must again calcu- 
late the work done in an infinitesimal expansion of the spherical shell. 
The system is in equilibrium when that work is zero. 

The forces in the present case are again derived from a potential, 
although this potential is not connected with the concentration of the 
oxygen alone. Since somewhere between r = 7, and r = 7, the force 
is zero, the potential has a minimum at that place. By an argument, 
similar to that which leads to equation (10), we find that the work 
done in an infinitesimal expansion of the shell vanishes when the po- 
tential of the effective force at r = 7, is equal to the potential at 
r = 7,.. But this means that the average gradient of the potential in 
the shell and hence the average force is zero. Approximately this 
means that the average concentration of oxygen in the shell is equal 
to that value c* , for which 6 = —B/A. Hence in order to determine 
the equilibrium configuration of the cellular shell we must determine 
the configuration, for which the average oxygen concentration c’ is 
equal to c*. 

The average concentration of oxygen c is determined again by 
considerations of material balance. Denoting by primes the constants 
that refer specifically to oxygen and remembering that for oxygen 
consumption q’ < 0, we have for material balance 

ae fev 5 
Ant.?D' ; (¢'2 — ¢'s) Last (25) 
fy 3 
which gives 


, , 


aq’ (2 — 11) 
cc, = ¢, — > 


3D"; 12? (26) 


Denoting by c’, the external concentration of oxygen and by D’, 
its external diffusion coefficient, we have outside of the shell (A, chap- 
teri, since 6 = 12) : 

Aar,?D', (C5 — C's) 


Ye 


4 
= a nag’ , (27) 
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which gives 


ia oe (28) 
The average oxygen concentration ¢’ is given by 
Wes an (29) 
From equations (26), (28) and (29) we have: 
Ge ee (30) 


Introducing now into equation (30) the expression (4) and writing 
the requirement of equilibrium c’ = c*, we obtain after rearrange- 
ments 
(2D'; + D’.) gar, — D'.q'av/r,3 —a = Lr? ; L = 6D',D'.(¢,— ©). 
(31) 
This equation determines 7, for every given a in terms of the other 
parameters. We have however the same difficulty as with equation 
(16). An exact expression for 7, in closed form is impossible. How- 
ever by the same argument as before we find the critical value a* of @ 
above which a cavity begins to form, by putting in equation (31) 
= Wa, and solving for a. This gives: 
6D',D’.(c¢'5 — c*) 
OF | eee 32 
ye: 
and 


6D',D'. (65 — c*) 
(2D, + D’.)q 


Putting Di; ~ D, ~ 107 cm? sec? ; gq’ = 10° gm cm® sec’, and 
ce) —c* wc, & 10+* gm cm we find 


(33) 


TAO Ose 


Ill 


The cases considered above are only two out of many other pos- 
sible ones. We may study in a similar fashion situations, in which q 
is a function of concentrations and obtain still different expressions. 
The interesting point of this investigation is that it shows that under 
rather general conditions formation of inner cavities will occur in cel- 
lular aggregates. It is the embryologist who should be most interested 
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in such studies, since formation of cavities of different sorts is an im- 
portant feature of embryonic development. 

It is of course highly unlikely, that all cases of formation of cavi- 
ties are due to the same mechanism. But such an assumption is not 
even necessary. In fact we already have found two different mecha- 
nisms which lead to formation of cavities and a number of other 
mechanisms may still be found. If we eventually wish to find the 
mechanism, responsible for the formation of one or another type of 
cavities, the natural procedure to follow is to study all theoretically 
possible cases and then by comparing the theory with observation de- 
cide which of the possible cases actually do occur. 

Consider for instance equations (18) and (33). While, as can be 
seen, they may roughly give the same value for the critical size at 
which the formation of cavity sets in, yet the dependence of this cri- 
tical size on various external factors is very different in the two cases. 
While equation (18) involves only the effective metabolism and the 
effective values of D and M,, equation (33) involves the correspond- 
ing values for oxygen specifically. If one could for instance observe 
the onset of the formation of a cavity in a morula, leading to a blas- 
tula, under varying conditions which affect the oxygen consumption, 
the different diffusion coefficients, etc., one might perhaps decide in 
favor either of (18) or (33). Thus, further mathematical studies in 
that direction will undoubtedly suggest new experimental problems 
to the embryologist. 

Speaking of the two cases discussed here, the case studied in sec- 
tion I seems to be more likely, particularly in the view of the fact 
that drugs inhibiting glycolysis do not seem to affect the formation 
of normal blastulae (Brachet, 1934). But one cannot entirely ex- 
clude the possibility of the other case. It is even possible that the two 
mechanisms may combine together. 

Inasmuch as in case of a negative effective metabolism no cell 
division should occur (A, chapter iii, MB, chapter x), and particular- 
ly no mistoses in superficial layers should take place one might think 
that the case studied in section II can have no embryological applica- 
tions. This conclusion is however not binding. In the outer layer the 
average q may be negative. But q will fluctuate from cell to cell, so 
that at every given moment some cells will have a positive q, and will 
divide. In the layer as a whole however the attracting forces will 
prevail. 

The following considerations may be of interest. In a single cell, 
considered as a sort of liquid drop, the formation of a cavity would 
require much stronger forces, because the inside of the cavity would 
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be filled only with the saturated vapor of the cell liquids with a very 
low pressure, whereas the outside surface would be subject to atmo- 
spheric pressure. But if we consider the cell as a sack formed by the 
membrane and filled with colloidal metabolising particles, suspended 
in water (MB, chapter viii; A, chapter iii), then formally each par- 
ticle will behave like a cell of an aggregate. These particles may well 
expand forming a sort of shell, with a cavity consisting of water free 
of particles and containing only nonmetabolising products. We would 
thus have a cell consisting of an outer metabolising layer and an in- 
ner, inert, watery layer. It may be interesting to speculate as to 
whether this may have some bearing on the theory of formation of 
vacuoles. 

The author is indebted to Mr. H. D. Landahl for checking the 
calculations. 
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A mechanism is discussed which has the property that two tem- 
porally separated stimuli presented over a common path may be 
discriminated as though they had been presented simultaneously over 
separate paths. 


In deriving some properties of a neural mechanism for discrimi- 
nation (Landahl, 1938), it was necessary to assume a simultaneous 
presentation of stimuli for discrimination to be possible. It is the 
purpose here to discuss a mechanism with such properties that two 
stimuli presented separately in time over a common path may be dis- 
criminated as though they had been presented simultaneously over 
separate but identical paths. The discriminatory response is made 
some time after the second stimulus is presented, but no response is 
made to the first stimulus alone. Such a mechanism, as is to be dis- 
cussed here, is not required, however, by a discriminatory mechanism 
of the type discussed by A. S. Householder (1940). 

Let us then investigate some properties of the mechanism repre- 
sented in Figure 1. We shall limit our discussion to the case in which 
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the mechanism is symmetric; that is, the members of any pair, as 
II and II, are identical. The following notation is adopted (Rash- 
evsky, 1938, ch. XXII). The symbol « denotes the excitatory factor 
produced at the efferent end of an excitatory fiber (solid line). The 
symbol 7 denotes the inhibitory factor produced at the efferent end 
of an inhibitory fiber (broken line). We are then assuming the fibers 
to be of the simple type so that each fiber produces either « or 7, but 
not both. Neither e nor 7 can be negative. The intensity of excitation 
E in a fiber excited by a stimulus S or by a net excitation e—j is 
assumed to be given by an expression 


E = B(S—h) 
or (1) 
E = B(e—j—h), 


where f is a constant of proportionality and h is the threshold of the 
fiber. Unless S > h we have E = 0. Arabic subscripts correspond 
with the Roman numerals. Where the Roman numeral denoting a 
fiber has a superscript, the symbols ¢« or 7 and EF associated with the 
fiber are primed; because of the assumption of symmetry, the other 
symbols which are fixed parameters of the fibers need not be primed. 
The superscript instead of a prime is used on the Roman numerals to 
avoid confusion with a convention used by Rashevsky (1938). 

The development of ¢ or 7 at the end of an excitatory fiber or at 
the end of an inhibitory is given by (Rashevsky, 1938, ch. XXII). 


de 

dt — AEH — QE ’ (2) 
adj ; 

ae BE — bj, (3) 


where A and a are parameters of an excitatory fiber and B and b are 
corresponding parameters of an inhibitory fiber. In Figure 1 the af- 
ferent end of each fiber at a synapse is represented by-a dot, and the 
direction of transmission at any synapse is readily determined. Ar- 
rows are also used to indicate the direction of transmission of nerve 
impulses in the fibers. 

We now introduce the following simplifying conditions to hold 
throughout the discussion 


A, = B;, a, = b;, he = hg. (4) 


; Let us assume that initially all e’s and j’s are zero. Then let a 
stimulus S > h be presented at the afferent end of fiber I. The rate 
of development of «’, and «, at s, and s, will then be given by an equa- 


ae 
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tion of the form of (2). Introducing (1) with the appropriate sub- 
scripts for fiber J and integrating the result, «’, and «, are obtained 
as functions of time, «’,(t) = e,(t) . We shall not need this expression 
but only its stationary value given by 


Ai 1 
faa = 8 ea (5) 


Since fibers IJ’, II, IIT’ and III are equally excited at all times and 
since j, = 7’. = js = 7's = 0 initially, we have e,(t) = e2.(t) = 7'3(t) 
= j,(t) because of conditions (4). Hence at s. and s; as well as at s, 
and s, the net excitation (« — 7) is in each case zero, &: — ge seca} 
— j, = 0. Under these conditions no excitation will be produced be- 
yond s, and s,. Now let the stationary state be reached to any degree 
of approximation so that after a time ¢ = 11 and until some external 


change is introduced, say at 12, we have 
A, Be 


go = & =—J3=—)s = (e, — ho) = 
ApBo A, pb : 
Se [AP s— iy —h | 
(6) 
Tail eth 


We assume that «, > h since otherwise «’. = & = J's = 9 and further 
discussion is meaningless. Likewise in subsequent equations where 
similar situations arise, similar assumptions are implied. 

Suppose now at a time t > 1 an amount of inhibitory factor 
Aj > 0 is produced at synapse s; due to some outside influence. The 
nature of this change will be discussed later. Since the net excitation 
acting on fiber IJJ is «, — 4 j instead of &, the inhibitory factor j; is 
decreased so that «’. — j; is increased. Since it was initially zero, it 
becomes positive. Let 


A> Ba 


2 


Aj>ha (7) 


so that at a finite time t; > t2 > 1, the increase in the net excitation, 
s', — j,, at synapse s, exceeds h, , or since e'. = &, 


fg > has te f= te, (8) 


where 1; is the time at which some change is again made, being the 
time at which S is to be withdrawn. The quantity t; — t may be ob- 
tained as a function of 47 but we shall not. need the value. Now, 
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after the time 13, 7’ > 0 because of (8) and thereafter Is and 7’, are 
interrelated. In the mean time it is clear from the mechanism that 
é’y , € and 7’; remain at their value given by (6), or 


Bearer! Ao Bs 


e Ss == 
2 2 73 de a, 


(Shs) — Ie 


(9) 
Ga = 005 tee Cite 


Since e, — 7's = 0, j4 remains zero. 

For the development of 7, and 7’, in the interval 7; > t > 13 we 
have the equations, (cf. Rashevsky, 1938, ch. XXIV, and Household- 
er, 1938) 


i = A; Bo (é1 Res A 7) SAG Bo qs eS Ae Ey ? (10) 
dj's ” 
“ap B+ Ba (62 — Ia) — Bs Bats — Oa (11) 


by substituting, in each case, equation (1) into (3), using the appro- 
priate subscripts, and introducing (4) and (9). The relations (10) 


FIGURE 2 


and (11) hold only if the respective E’s are positive, or rewriting, we 
have the conditions 
js <e:—h,— Aj, orelsej, =0, (12) 
Js < &—h,, orelsej’,=0. (13) 


At the time t= 7, , 7’, = 0 so that condition (12) is satisfied. Because 
of (8), condition (18) is also satisfied initially. 
We now wish to determine sufficient conditions which insure that 
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after a sufficient lapse of time (but still within the interval 1; >1 
> 13), 7; has a value zero while j’, has some positive value. Setting 
the derivatives in equations (10) and (11) equal to zero and solving 
for j’, in each case, we determine two straight lines: (Figure 2) 


a, é Com. 
a: de ea ae AL (14) 
} B By, Bs . 
LE ersA Bia ae (€2 — Is) aeseet y (15) 
on which respectively dj,;/dt = 0 and dj’./dt = 0. 
Now let 
fa te fo O.. (16) 


Then because (8) implies that «. > h,, the lines Ls and L’', have posi- 
tive intercepts and negative slopes. If L’, has a larger intercept on 
the j,-axis than L, and if L’, has a larger intercept on the j;-axis than 
L, , that is, if 


B, 4 
é:—h,— AJ =! ue (& — Nx) ? (17) 
c Ae 
pete Aig (ela) 3 (18) 
which may be written, because of (9), 
As Be 4 B, 4 hg 
ee tee) (i ee ee (19) 
ops bs 
. phy 
Lara hs 


then L’, and L, will not intersect in the positive quadrant. The posi- 
tive quadrant (j,; and j’, must be positive to be meaningful) is then 
divided into three regions; J, the region below the line L; ; IJ, the 
region between L, and L’; ; and III, the region above L’, (Figure 2). 
The points on the axes are included in the respective regions but not 
the points on the lines L; and L’, which separate the region. 

In region I, 


dj’, djs f 
i >0 and St OPS 


in region II 


dj's > 0 and djs 


at pet 
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and in region JJI 


dj's dj; 
at < 0 and Ti cal Vie 

However these inequalities are subject to conditions (12) and (18). 
These conditions are satisfied in the region bounded by the broken 
lines of Figure 2. A point to the right of the vertical broken line 
fails to satisfy condition (13) and hence j’, must decrease to zero if 
somehow 7; were maintained. But dj7;/dt = 0 and hence 7; decreases 
until the configuration point enters the rectangular region required 
by (12) and (13). A point above the horizontal broken line fails to 
satisfy condition (12) and hence 7; must decrease to zero. However 
if 7, = 0, from (11) we still have dj’,/dt > 0, and the configuration 
point moves to the stable value given by the larger ordinate inter- 
cept, L’, with the 7’,-axis. Now the only possible stable equilibrium 
configurations are the four intercepts, since clearly the point (0, 0) 
represents an unstable equilibrium, and since from the nature of equa- 
tions (10) and (11) and since 7, > 0 and 7, > 0, neither 7; nor 7’, can 
increase without bound. From an inspection of Figure 2 it is clear that 
the intercepts of L; and L’, on the 7;-axis are unstable points of equi- 
librium as is also the intercept of L; on the 7’,-axis. There remains 
the intercept of L’, on 7’,. Because of condition (12) 7; = 0 and re- 
mains so. From equation (11) with 7; = 0, it is clear that a displace- 
ment of 7’, away from the intercept value gives such a sign to dj’,/dt 
that the initial intercept value tends to be restored. Let this equilib- 
rium be reached to any degree of approximation at the time t = 1, 
™, > 13 ; we may then write 


*f era 
Pr tie A, 


(5 = h,) ’ 


(21) 
js =0,forz; >t>u. 


In the above argument we have imposed four conditions: (4), 
(7), (16), and (19), since (7) implies (8) and is identical to (20). 
Now condition (4), though a very restricting condition, simply re- 
lates certain parameters in two pairs of fibers. Conditions (7) and 
(19) impose a minimum value for Aj while condition (16) imposes 
a maximum value for 47, a value which increases with S. If the 
right hand member of the inequality (20) is greater than the right 
hand member of (19) then (7), (16) and (19) may be written 


Ash, 


ee Peg Set Ang ee 
€ > tere Rak (22) 
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and hence, since €, 1S a monotonically increasing function of S, (22) 
imposes a minimum for the value of S greater than the threshold 


h,, or because of (5) 
A, A, he 

A, Az Bi Be 

But if the right hand member of inequality (20) is less than that of 
(19), then (7), (16) and (19) may be written, instead of (22), as 


: Az Be Bs Bs By, Baha 
eh) (24) 


MCS 
S> hy pag + (23) 


or since the first member must exceed the last member, we again 
have, after simplifying, 


a Ne Ay Az hy 
Sc ARE Ay Ang, pas 


Hence if we assume conditions (4) and (25), then if 47 lies within 
the range given by (22) or (24), and such a range exists, the result- 
ing stable values of j,; and j’, are given by (21). 

Now let the stimulus S be maintained after equilibrium has been 
reached and long enough so that the set of self-exciting neural cir- 
cuits, represented in the figure by C’, which have sufficiently low 
thresholds, are brought into a state of permanent excitation (Lan- 
dahl and Householder 1939, Rashevsky 1938, ch. XXIV). For sim- 
plicity, let us assume that the parameters are such that for. any cir- 
cuit the threshold value for permanent excitation to be maintained is 
approximately equal to the threshold h, of the fiber. Also we assume 
that all the parameters of the circuits represented by C’ are the same 
except the thresholds, which are distributed according to some dis- 
tribution function. Now all those circuits of C’ are excited for which 
he < @(S). In the meantime, no fibers of C are excited.. If we as- 
sume that the distribution of thresholds he is such that in each range 
A hz there are an equal number of circuits, then, as the circuits are 
otherwise identical, the net excitation ¢'c acting on the afferent ends 
of fibers VII‘, VIII‘ and IX' is proportional to the number of cir- 
cuits excited or A 


(25) 


eo =De,(S), 8 “7 (26) 


where D is a constant. 

-- -Now any time after the circuits have become self-exciting the 
stimulus S , as well as 4 j , may be withdrawn. We then. consider the 
effect beyond C’ at t > 15. If we let hs < h,, Bibs > Abc, and 
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Bzfs/bs > Azb:/@,, then at all times j; > «, so that fiber X? is not 
excited. Let the stationary state be reached to any degree of accur- 
acy after t > t. > 1s, So that we have the following relations, 


A, ue 
&; = P (€'¢ =—= h-) (27) 
a, 
B 
7 SE Cp ey (28) 
bs 
Bs Bo 
diners a (e'¢-— he), (29) 
Jige>- Oe (30) 
for 
Tz > t => T6 o) 


where 1, is the time of a new external change, being the time of pre- 
sentation of a new stimulus S. The value of 7’; in (30) is not given 
for reasons given in the next paragraph. Now in the interval 1, > f 
> 7,, all e’s and j’s are zero except eo, &,, 9's .J9,andj,;. But after 
S is presented, t > 1,, €, and «’, develop to their stable values given 
by 

evans (Sey (31) 


1 


Similarly, & , e', and j; develop to their stable value given by (since 
Aj =0) 


, 


t2—— 16 oF) a 


A, ze By 


As Cy, 


(S—m) —h], (32) 


with «’, — 7, = 0 at s. and s, throughout the development. But 7’; 
cannot develop at the same time, because of 7’;. Hence the net exci- 
tation at s; is positive and fiber IV will be excited if 7’; is sufficiently 
large. We then have j’, and j, interrelated, during the stationary 
state of e’, etc., in the same manner as j,; and 7’, were in the interval 
ts; > t > t3 but with 7’; playing the role of 47. If 7’; satisfies the con- 
ditions satisfied by 47, then 7’; tends to its stable value, zero. Hence, 
at time z,, after sufficient time has lapsed, the circuits C, whose 
thresholds are less than e, , are excited, and we have 


& = De,(S), fort > 73. (33) 


Any time after ts, stimulus S may be withdrawn. 
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If fibers V and V: were such that for them we may use equation 
(1), the above argument is complicated and not strictly valid since 
yj; is a function of S. In order to avoid such a difficulty let us assume 
that these fibers are such that H;(S) is given by a simple exponential 
[Rashevsky, 1938, ch XXII, eq. (5)] and the parameters are such 
that for the S and S used, 7; is greater than 0 and essentially inde- 
pendent of S. 

Consider now the situation after 1), when é¢ has acted sufficient- 


ly long for é;, jg and 7, to reach their stable values with «, and j, be- 
ing given by 


ee A: be 
bia (34) 
a, | 
comet gt 8g if 
jo = —— (6 — he), (35) 
9 
t = Tg « 
Now «’, and j’, are unchanged and hence also 
A 
yeaa Cee ay (36) 
Az : 
B | 
jo= zi (eo — he) , (37) 
9 
; t > Tg. 


However j’, is no longer given by (28). If ec — jac Nathatis it 


(38) 


j,=0, te te (39) 


Since S ~ § in a discrimination situation, it is sufficient, for (38) to 
hold, that B.f./by > 1 if hy < Dohs/Bof.. But then also, at least if 
Bzf/b, >> 1', 


by — bo h 
Eg — B. i Ec > hy Ta 2 (40) 
and hence 
Ec et 4% mS hs ’ 
and 


je=0, td. (41) 
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Therefore at t = 1), fibers X! and X are respectively excited by 
the quantities «',(S) and «,(S), of equations (36) and (34). How- 
ever, during the interval 1, > t > 1, fibers X1 and X are not neces- 
sarily excited equally. We may resolve this difficulty in two ways. 
We may assume that e; and 7’; develop much more slowly than j, or 
é> so that e- and 7’, may be assumed suddenly established at t = 7. 
Then if the parameters of VIJ and VIII are related by conditions 
analogous to (4), the development of the net excitation at s, and s; 
will be simultaneous. On the other hand we may disregard any tem- 
porary differences at s; and s; and assume that the next higher order 
fibers, X' and X are of such a type that a,. is much larger than @,, 
b,, or b, and hence the development at the net excitations at s; and 8; 
are essentially simultaneous so far as X! and X are concerned. 

Suppose that X! and X lead to a discriminatory center of the 
type previously described (Landahl, 1938). The value 


; ca A, B; ae A, B; Az Bs A, 6, D 
er — eer (e'¢ — &¢) =——_____—_ 
A; Az Ay Ay 


(S— 8) #3 


is analogous to «, é. of equations (5) and (9) of the previous de- 
velopment. This latter value was proportional to the logarithm of the 
ratio of the two simultaneously presented stimuli, but since they are 
necessarily very nearly equal, the logarithm of the ratio is approxi- 
mately proportional to the difference. If instead of the first equation 
of (1) we had used, as in the previous case, a logarithmic relation be- 
tween E' and S, then in place of (42) we would have e’, — e,; propor- 
tional to log (S/S). 

Under the conditions imposed, the mechanism of figure 1 has the 
property desired. As far as fibers X' and X are concerned, or any re- 
sponses that they may lead to, it makes no difference in the response 
whether two stimuli S and S are presented in succession at fiber J, 
each being sufficiently large and being presented for a sufficiently long 
time, or two other stimuli «’;(S) and «,(S), which are linear func- 
tions of S and S, are presented simultaneously at s, and s; to X‘ and 
X directly. 

In order that the mechanism of figure 1 may discriminate re- 
peatedly between pairs of stimuli, suppose also that the response to 
either S or S leads to an inhibition at synapses s,, and s, sufficiently 
large, and for a sufficiently long duration so that «’, and é¢ subside 
to zero. Then all e and 7 not already zero will decrease to zero. A 


second pair of stimuli may then be presented and the process re- 
peated. 


— 
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A large number of conditions have been imposed on the para- 
meters of the mechanism. Not all have been necessary. To eliminate 
those which are not necessary would require a more detailed analysis 
augmented by a comparison with experimental data. Such an assump- 
tion as perfect symmetry of the mechanism is certainly not necessary 
since temporal and spatial biases are generally present. However, by 
appropriate technique, they may often be made neglible if desired. 
Also it was assumed that the distribution curve of hy was a constant. 
Now almost any distribution function having one maximum which 
is not too narrow could have been assumed without changing the 
results essentially. Also it was assumed that the threshold for per- 
manent excitation for each circuit of C’ and C was approximately 
equal to the threshold of the fiber. Such would generally not be the 
case. If not, the second stimulus would be favored so far as this one 
factor is concerned. 

It was assumed that 47 was somehow applied at synapse sa. 
Without much change in the results, 47 may have been applied at 
$2, or at Ss, or s , though in the latter cases, the responses would have 
been reversed. But this would be satisfactory. Instead of 47 > 0,a 
quantity A ¢ at s,, S), Sc, or 8a could have been used. One might have 
assumed a slight asymmetry in the mechanism such that leg lie OY 
h’, ~h,. We consider Aj (or its equivalent) as some random change 
and this change must exceed a certain value before the initial system, 
in which «’, = 7;, passes over into a much more stable system in 
which e’, is unchanged and in which 7; = 0. Thus 47 measures the 
half width of the stability of the initial configuration «2 = es = Js 
= j’, of equation (6). 

The author is indebted to Dr. A. S. Householder for discussing 
and checking this paper. 

This investigation has been aided in part by a grant from the 
Rockefeller Foundation to the University of Chicago. 
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Laminar motion of two viscous incompressible fluids through 
each other is treated for two cases: flow along the axis of a circu- 
lar cylinder, and flow between parallel flat plates. Motion of either 
fluid entails that of the other. Regarding one fluid as a solvent, the 
other as a solute, and supposing the system to have ends imperme- 
able to the former, it is found that the solvent streams with the so- 
lute down the center of the system, to return in the opposite direc- 
tion out nearer the walls. Thus diffusion of a dissolved substance 
through a region in which the solvent is confined produces continual 
streaming in the latter. 


I 

Consider two fluids moving through each other under the follow- 
ing conditions: (1) The fluids have constant densities. (II) The 
motion is slow enough that inertial forces may be neglected. (III) 
The steady state has been sensibly reached so that time variation 
terms may be neglected. (IV) Long range forces are negligible, and 
the short range forces between the fluids may be represented by a 
drag force proportional to their relative velocity. 

Letting P denote the mean pressure of a fluid, H its coefficient of 
viscosity, and (U, V, W) its velocity, with corresponding small let- 
ters for the other fluid, the equations of motion are (Young, 1938) 


P, = B(u—U) + HV2U 
= B(v—V) + HV°V 
P, = B(w—W) + HV°W 
Pr = B(U—Uu) + Vu 
py = B(V—v) + V0 


p. = B(W—w) + Vw (1) 
OU OV OW 

Spulaayes page 

Cu ov ow 

Gala hi Gh TA Uae eS 9 | 

on ee Oz 
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Herein P, = 0P/dx etc., and B is a constant. Differentiating and add- 
ing gives 
V?P=V?p=0, (2) 
so that the pressures are harmonic functions. By addition in (1) 
P,+ pe = HV?U + Vu (3) 
and similar equations; while by differentiation arise such equations as 
V2U = Viu—F Vtu , (4) 
from which Vp, = 0 has been dropped. Substitution of (4) into (3) 
gives the system of equations 
Viu+ aV*Uu = Ge 
Viv + aV?v = Dy 
Viw + aV?w = gz 


(5) 
ou ov ow 
aig ak oe a 
where 
(H+ 7)B 
ae Hy 
(6) 
1 2B) s 
—= To : 


The corresponding solvent velocities are found from (1) to be 


yt a 
U=4u BVUT sR 
V=v—3yv+4 (7) 


==, Sagar 2 D: 
W=w es 


which may be readily verified to satisfy (1) for any u,v,w satisfying 
(5). Thus the most general solution of (1) is given by the general 
solution of (5) supplemented by (7). 

In what follows attention will be restricted to systems in which 
the fluids move only in one direction, in which case the above rela- 
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tions are considerably simplified. If the x axis is taken along the di- 
rection of flow, it follows that v = w= V = W=0. From (1) it is 
then seen that P, = P, = py = pz = 0, so that the pressures are func- 
tions of x only. But (2) then shows that they must be of the form 


P=kz+M 
(8) 
p=kxe+m. 
The system (5) then reduces to 
Vtu+aV*u=ob 
(9) 
ou _ 
ee as 
where 
» (K+ k)B | 
bie serra (10) 
while (7) becomes 
ee rege ee 
U=u—z Vit =e (11) 


It may be noted that restriction to one-dimensional or laminar 
’ flow renders the initial assumption (II) superfluous. For the terms 
thereby neglected in (1) are of the form 


(12) 


and vanish identically because of (9) and the condition v = w = 0; 
the same is true, of course, for the other fluid as well. This would not 
be true for a compressible fluid, since it would experience finite accel- 
erations as its pressure changes along the path of motion. 


II 
Consider first the case where the flow is along the axis of a cir- 
cular cylinder, assuming symmetry about that axis. Then u = u(r) 
and 


du idu id, du 
\/ = ——+ -— = -=— ——') ¢ 13 
= aT pdr. rar Orie ey 


and similarly for U. Equation (9) then becomes 


ad, dV*u 


d , du 
ate ke ee lpn) ee 14 
4p (r a )-+as (7) br . (14) 
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The integral is apparent and upon division by r is 


ad yh 18) A, 

ee 2 Se SS — 15 

RAE AOA et ip Sse f (15) | 
where A, is an integration constant which must be equal to zero to 
keep u regular at r = 0. A second integration gives 


Viutau=F99+ Ay ; (16) 
or 
du idu b. 
Ga eaceie ) get tere — 


The homogenous equation obtained by setting the right side equal 
to zero is a zero order Bessel equation with solution 


AzJ.(Va@r) + A, Yo(VGr), (18) 
where J, is a Bessel function of the first kind and Y, is a Neuman 
Bessel function of the second kind. Since Y,(0) = 0, A, mustine 


equal to zero. 
It may be verified by substitution that 


Cee As b 
a+ -s (19) 
is a special solution of (17) ; hence the general solution is 
oie b 
UE AGT A NIRT Lee Sale (20) 
where 
A, b 
A= we (21) 
From (11), along with (16), (20), and (21) it is found that 
an A b k = »b 
U = (1 fe ae es ee 
( rep) Aes 6G CT et es areata ane (22) 


Both fluids have zero velocity at the surface of the cylinder r = R. 
Writing u= U=0forr=R gives two equations for determining the 
constants of integration. They come out to be 
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rcp 
pees ee a 
and (Va R) 
(23) 


and with these values (20) and (22) become 


u = E(r?— R?) +H A[CJ.(Var) —1], 
(24) 
U = E(r? — R?) —7 A[CJ,. (Var) —1] 
where 
_ K+k 
4(,+H) 
kH — Kn 


i\ 3 2S 
By +H)? (25) 


a 1 
Jo(Va@R)~ 


Since a is negative the Bessel functions appearing in the above 
have imaginary arguments, and hence do not vanish for any value 
of r exceptr=0. 

The functions 7? and J,(,\/ar) are linearly independent; hence 
the flow rates in (24) can vanish identically only if EH = A= 0. 
From (25) it is seen that H = 0 requires K + k= 0, and then 
A = 0 implies K = k = 0 . Hence if there is a pressure gradient in 
either fluid then both will be in motion. Regarding one fluid as a sol- 
vent and the other as a solute, this shows that any diffusion of the 
solute is accompanied by convection currents in the solvent. 

The total flow thru a cross-section of the cylinder is 


R ‘ Die 
ya2af pee — 229A [oR VE_- |e 
; VaR) 2 


(26) 


R 4 Ea 
J =2n | He peace +2aHAl oliver) = tae 
2 VaR) 2 


0 


Supposing & given, then (26) may be solved for that value of K 
which makes J = 0. Thus if we imagine a system closed at 2, and 
x. by membranes impervious to the solvent (capital letters) but al- 
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lowing penetration of the solute (small letters), and if the solute is 
under different pressures p, and p, at the two ends, then there will be 
a net flow of solute thru the system but J will be zero, K being there- 
by determined. Since the solvent cannot be everywhere at rest it 
therefore will flow in one direction in some parts of the cross-section 
and in the opposite direction elsewhere. There will thus be a cir- 
cuital streaming of the solvent; as will be seen later it flows with the 
solute down the center of the cylinder and returns in the opposite. 
direction out nearer the walls. It has to turn at the ends and there the 
flow will not be confined to the x direction, so that the above treat- 
ment does not apply. However if the system is sufficiently long the 
above considerations hold, except near the ends. 

Thus we have an example wherein a substance diffusing thru a 
region in which the solvent is confined sets up continual streaming 
in the latter. 

In typical diffusion problems the constants are such that the 
above expressions may be considerably simplified. The bracketed ex- 
pression in (24) is equal to —1 at the center of the cylinder and rises 
monotonically to zero at the walls. If \/—aR (which is a dimen- 
sionless quantity) is large enough the bracket will be sensibly —1 
everywhere except very near the walls where it increases rapidly to 
zero. Except in this transition layer itself the Bessel function term 
may be ignored. 

For x fairly large J, (iz) = iJ, (ix), so that the bracket term in 
the flow rate (26) is approximately 

if 1 


Tak ase (27)4 
Calculation for typical biological cases (Young, 1938) shows that 


V—a is of the order of 10° to 10° cm~, so that if R is say 10+ cm. or 


more the first term in (27) is negligible compared with the second. 
Thus, very nearly, 


aH 


J =— 5 R* + a7 AR? 


_ (28) 
} = — "= Rta HAR? 


Setting J = 0 for the case of a cylinder closed to the solvent at its 
ends gives 


i A (29) 
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and in this case (28) gives 
j =—a(y+ H)AR?. (30) 
From (29), (25) and (6) 


8 7 a 8 
K+ aye | =—*| 1— ae]: (31) 


and with (—a) R? of the order of 108 or more as indicated above it is 
seen that K may be taken as —k in A to give very nearly 


k 


A = ———_~. . oa 
BC+ 1) So 
The total solute flow (30) is then approximately 
: k 
4 = ——a2ak’. (33) 


B 


But (34) is exactly the total rate at which the solute would be trans- 
ported by simple Fick diffusion in a stationary solvent, so that the 
solvent streaming has very little effect upon the total solute trans- 
port down the cylinder. 

For x large Jo(ix) = e*/\/222, and in the velocity expressions 
(24) 


CIo(Var) = , ee (34) 


which for \/—a large is entirely negligible except very near the outer 
wall of the cylinder. Thus everywhere except in the wall layer (24) 
and (32) give 


k 4 2r?— R? 


Meas H+ Te 
(35) 
ke 2yr? — (H+ 2) R? _ ie 


Be: (H + 7) R? Bs 


The term — k/B in wu is the velocity of simple Fick diffusion in the 
stationary solvent. 

Down the center of the cylinder, out to r = R/v2 the solvent flows 
in the same direction as the solute; outside this central stream the 
solvent flows back in the opposite direction. The cross-sectional areas 
of the two streams are equal, and the total flow in either is given by 


 —— 
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R/ 42 ken eae 
— See ge tlla ff7, 36 
J 2a [ rUdr eee (36) 


Dividing by the cross-section area of the stream gives the average ve- | 
locity in it 


pal 
(cen ee (37) 
B2H+y7 
From (33) the average solute velocity is 
k 
ee poe 38 
Uw rao (38) 
so that the ratio of average “streaming” to “diffusion” speed is 
Ue 
w 2H+y 


From (85) or, more exactly, from (29), (32), and (24), it is 
seen that the velocities remain finite even when R increases indef- 
initely. 


III 
Similar results are obtained for the case of laminar flow between 
parallel flat plates. Let the y axis be perpendicular to the plates 
which are at y = 0 and y = L respectively, the x axis being along 
the direction of flow as before, and assume the flow independent of 
z. Then u = u(y) and (9) becomes 


d‘u d?u 
—+a—=b. 40 
Writing Z = 2 this i 
g => dy 18 1S 
CZ 
dy? +azZ=b, (41) 
of which the general solution is 
du e _ b 
a a evevt A, aie Sayre : (42) 


Integrating twice gives 


u— A, v-ay A, -V=a y b 2 A A.: 
ar “TParery aimeree ae A ie sy + ra 9 (43) 
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and with (11) and (42) 


U = — A= 44) oA, 4 
(44) 
nb 
Pa 
The boundary conditions u = U = 0 at y = 0 and at y = L then 
yield, denoting 1 = ~—aL, 


b k 
anor + Ad + Arta 


eee re al 
(45) 
bL 
Drei a gb 
4 a7 a 


Introducing these values into (43) and (44) and making the change 
of variable ~ 


> (46) 


so that Y is measured from midway between the plates, there is final- 
ly obtained 


U= F(y?—=) — 7 A[G cosh(V—aY) —1] 


(47) 
2 — 
w= FY? a) + HA[G cosh (./—aY) —1]. 
Here A is given by (25), while 
hies2h 
(48) 
e!/2 pores e-!/2 
e et — et 


As in the cylindrical system U or u can vanish identically only if 
both K and k do, and hence diffusion necessarily produces convection 
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currents in the solvent. The total flows through a cross-section of 
unit length in the z direction are given by 


L/2 3 Zi , l 
ippes vay = — Aa] —— G sinh — — L 
-L/2 6 ~ 
(49) 
L/2 3 l 
= uy =~" + AH[ 2G sinh — 1 |: 
ae 6 Ve 2 


Setting J = 0 and thus determining K gives the case of a system 
closed at its ends (a, and x.) to the solvent, just as before. 


For \/—a L large the bracket in (49) reduces to 


2 (50) 


in which the first term is negligible. Then, very nearly, (49) becomes 


ae = L' +7 AL 
(51) 
= — L? —HAL. 
Putting J = 0 gives 
6 

| eas 73 A’, (52) 

or, with (6), (48) and (25), 

12 n 12 
K =a 1 aT Sein | beset ee catemeeen| _—_— 

=e | x] 1 —an| (53) 


For (—a)L? large, A again reduces to the value given in (32); and 
then from (51) and (52) 


Ee k 
Saas DL, (54) 
which is the simple Fick transport through the stationary solvent. 
For | large the quantity G of (48) is approximately 2e-'/2, and 
in (47) 
G cosh (/—a Y) = e-“ahre» 5 (55) 


which is negligible except near the walls. Except very near the walls 
the velocities are given by . 
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yn 12Y?—L? 
n 21? 


(56) 


The solvent velocity changes sign at Y? = = L? , so that the 


central stream is wider and slower than the returning stream. The 
total flow in either direction is 


Jt=2 OU aY #01926 2 SL, (57) 
: oer 1: ee eas 
the mean speed in the central stream being 
feo n 
[Ji Ceo ha het, Ee lene Ale een 58 
U Fae: ee (58) 
and in the outer streams 
k n 
pe ng | ee 59 
U = 0.457 Peay (59) 


The mean solute velocity is again given by (88) and as before the 
ratio U*/u* is of the order of (89). 

As before, the velocities remain everywhere finite even if L be- 
comes infinite. 

This investigation has been aided in part by a grant from The 
Rockefeller Foundation to the University of Chicago. 
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